arXiv: 1507.06799vl [cond-mat.mtrl-sci] 24Jul2015 


Theory of perpendicular magnetocrystalline anisotropy in Fe/MgO (001) 

Dorj Odkhuu 1,2 , Won Seok Yun 1,3 , S. H. Rhim 1 ^ and Soon-Cheol Hongfj] 

1 Department of Physics and Energy Harvest Storage Research Center, 

University of Ulsan, Ulsan, 680-749, Republic of Korea 

2 Department of Physics, Incheon National University, Incheon, 406-772, Republic of Korea 

3 Department of Emerging Materials Science, DGIST, Daegu, 711-873, Republic of Korea 

(Dated: July 27, 2015) 

The origin of large perpendicular magneto-crystalline anisotropy (PMCA) in Fe/MgO (001) is revealed by 
comparing Fe layers with and without the MgO. Although Fe-0 p-d hybridization is weakly present, it cannot 
be the main origin of the large PMCA as claimed in previous study. Instead, perfect epitaxy of Fe on the MgO 
is more important to achieve such large PMCA. As an evidence, we show that the surface layer in a clean free¬ 
standing Fe (001) dominantly contributes to E M ca> while in the Fe/MgO, those by the surface and the interface 
Fe layers contribute almost equally. The presence of MgO does not change positive contribution from (xz\£z\yz), 
wherease it reduces negative contribution from (z 2 \£x\yz) and (xy\£x\xz,yz). 


Exploration for magnetic materials with future applications 
dates back more than two decades, which includes giant mag¬ 
netoresistance (GMR) and many applications in spintronics 
such as magneto-resistive random-access-memory (MRAM), 
magnetic sensors, and novel programmable logic devices LU1. 
Among those, perpendicular magnetocrystalline anisotropy 
(PMCA) has attracted greatly as materials with large PMCA 
can offer more opportunities to realize magnetic devices. It 
can provide an ideal tool to realize spin transfer torque ex¬ 
cluding external fields. Also it offers large bit density in prac¬ 
tical applications. Magnetic tunnel junctions (MTJs) using 
MgO as an electrode, in particular, have drawn huge atten¬ 
tion owing to enormous magnetoresistance reaching as high 
as 180%|2i]. Furthermore, large PMCA in Fe/MgO has shown 
great promise with tremendous performance |3 h5|]. 

There have been numerous subsequent works ever since in¬ 
cluding Fe/MgO and FeCoB|MgO. In theoretical sides, sev¬ 
eral ways have been proposed to enhance PMCA in Fe/MgO: 
By applying an external electric field®®, adding heavy tran¬ 
sition metal layer on top of Fe®, and so forth. While it is 
still on-going endeavor to enhance PMCA, it looks quite indis¬ 
pensable to identify the physics origin of PMCA in Fe/MgO 
without aforementioned effects. Although it has been at¬ 
tributed to the p-d hybridization between Fe and O atoms at 
the interface!9, l|] that is responsible for large PMCA, we 
argue that the hybridizationcannot be the main driving force. 
If the hybridization is indeed the main source of PMCA, why 
PMCA is absent in Fe grown on A10 x |Qjl[l2] despite the pres¬ 
ence of apparent hybridization between Fe and O? Further¬ 
more, if the hybridization is really the main origin of PMCA, 
then it should be small when Fe/MgO is underoxidized due to 
the reduction of hybridization ||. 

In this paper, we show that the most dominant contribu¬ 
tion to PMCA is not the p-d hybridization despite its weak 
presence in Fe/MgO. Instead, perfect epitaxy of the interface 
would be the key, which is an indicator of high level of sample 
fabrication. To support our idea, the contribution from the sur¬ 
face and the interface Fe in the presence of the MgO is com¬ 
pared, which is almost the same in magnitude. This clearly 
counter-argues that the hybridization plays the main role of 


PMCA, whose detailed analysis of the electronic structure is 
provided. 

Density functional calculations are performed using the 
highly precise full-potential linearized augmented plane 
wave (FFAPW) methodil3l] and Vienna Ab-initio Simula¬ 
tion Package (VASP)|[l4, Hfl. Generalized gradient approx¬ 
imation (GGA) by Perdew, Burke, and Ernzerhof (PBE) 
parametrization[ ITS is employed for the exchange-correlation 
potential. In VASP calculations with projector augmented 
basis sets lfnll . full atomic relaxations have been carried out 
within force criteria 0.001 eV/A, where cutoff energy 400 eV 
for wave function expansion and 16xl6xlk mesh in the irre¬ 
ducible Brillouin zone wedge are chosen. In FFAPW calcula¬ 
tions, cutoffs for wave function and potential representations 
are 16 and 256 Ry, respectively. Charge densities and poten¬ 
tial inside muffin-tin (MT) spheres were expanded with lattice 
harmonics £ < 8 with MT radii of 2.1, 1.4, and 1.8 a.u. for Fe, 
O, and Mg atoms, respectively. For k point summation, we 
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FIG. 1. (color online) (a) MgO-free Fe layers and Fe/MgO with sub¬ 
strate with (b) 1-ML, (c) 3-ML, and (d) 4-ML Fe. Spheres with yel¬ 
low, red, and green color denote Fe, Mg, and O atoms, respectively. 
S and I denote atoms at the surface and the interface, S-l that of 
subsurface. 
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TABLE I. Magnetic moments (in ps) of Fe atoms in the clean Fe and Fe/MgO, where Fe(S), Fe(S-l), Fe(C), and Fe(I) stand for Fe atom at 
surface, subsurface, center, and interface, respectively. Without MgO, Fe(S-l) and Fe(C) are not well defined for 1-, 2-ML, so is Fe(S-l) for 
3-ML. 


# of Fe layers 

clean Fe 

Fe/MgO 


FefS) 

Fe(S-l) 

Fe(C) 

Fe(S) 

Fe(S-l) 

FefI) 

O(I) 

1 

3.14 

- 

- 

- 

- 

3.09 

0.03 

2 

2.82 

- 

- 

2.85 

- 

2.62 

0.02 

3 

2.96 

2.39 

- 

2.96 

2.45 

2.83 

0.02 

4 

2.94 

2.45 

- 

2.92 

2.44 

2.72 

0.02 

5 

2.96 

2.43 

2.62 

2.95 

2.44 

2.80 

0.02 

6 

2.95 

2.44 

2.56 

2.95 

2.44 

2.79 

0.02 

7 

2.95 

2.43 

2.42 

2.95 

2.44 

2.79 

0.02 


used 24x24x 1 mesh in the irreducible Brillouin zone wedge. 
A self-consistent criteria of 1.0 x 10 -5 e/(a.u.) 3 was imposed, 
where convergence with respect to the numbers of the basis 
functions and k points was also seriously checked. MCA en¬ 
ergies ( Emca ) are obtained using FLAPW. For the calculation 
of Emca , torque method[18] was employed to reduce compu¬ 
tational costs, whose validity and accuracy have been success¬ 
fully proved in conventional FM materials and others [Hi¬ 
ll 

The model geometry in our study is depicted in Fig. |TJ (a) 
MgO-free Fe-layers and (b-d) Fe-layers on the MgO substrate, 
where we have considered number of Fe layers from 1 to 7 and 
five layers of MgO. The experimental lattice constant of MgO 
(4.214 A) was adopted for the in-plane lattices. Magnetic mo¬ 
ments of Fe atoms are listed in Table [fl where Fe atoms at the 
surface, subsurface, (one layer beneath the topmost surface), 
and the interface with MgO, are denoted by Fe(S), Fe(S-l), 
and Fe(I), respectively, as in Fig. Q] The center layer in the 
absence of MgO is labelled as Fe(C). Since Fe(S) and Fe(I) 
are identical 1-ML Fe/MgO, only Fe(I) is shown. In the clean 
Fe layers, i.e. MgO-free Fe layers, moments of Fe(S) are 
enhanced with respect to the bulk giving maximum value of 
3.14 iab for 1-ML Fe. Moments of Fe(S-l) are comparable 
or smaller than those of Fe(C). Fe(C) moments are slightly 
larger than bulk Fe, which is due to enlarged lattice constants. 
The presence of MgO also enhances moments of Fe(S) up to 
2.96 iab . In all cases, Fe(S) have the largest moments mostly 
around 2.95 /!#. Interestingly, Fe(I) moments are smaller than 
but comparable to those of Fe(S). This feature has been al¬ 
ready addressed in previous work^b]. 

Emca as function of the number of Fe layers is plotted in 
Fig. [2| with and without MgO. The presence of MgO systemat¬ 
ically increases total Emca by ^20% with respect to the free¬ 
standing Fe-layer except 2-ML Fe. As seen clearly, the satu¬ 
ration behavior is evident as the number of Fe layer increases. 
The singular behavior of 2-ML Fe will be discussed later. To 
reveal the role of Fe 3t/-0 2 p hybridization, Emca is plotted as 
function of the Fe-0 distance in Fig. Ob) for the case of 6-ML, 
where the dotted line denotes Emca of 6-ML Fe without MgO. 
Emca decreases as the Fe-0 distance increases. When the dis¬ 
tance exceeds 3.0 A, MCA energy becomes equal to that of 
MgO-free Fe layers. Definitely, this implies that hybridiza¬ 


tion affects Emca • However, in forthcoming discussion, we 
will show that the hybridization, though not completely ig- 
norable, is not the main contribution to PMCA as claimed. 

The atomically decomposed Emca is presented as function 
of number of Fe layer in Fig. (2c) and (d) for the clean-Fe and 
Fe/MgO, respectively. In the clean-Fe, the surface contribu¬ 
tion is largest, while those from the subsurface and the center 
layers are much smaller. Contributions from the subsurface 
are even negative. On the other hand, with MgO the interface 
layers in contact with the MgO contribute almost equally as 
the surface layer when the number of Fe layers exceeds three. 
Contributions from layers other than interface and surface are 
negligibly small. [See Supplementary Information]. If the p-d 
hybridization is really the main source of the large PMCA ob¬ 
served in the Fe/MgO, then why the surface Fe layer with no 
hybridization contributes almost equally as Fe(I)? This ques¬ 
tion will be explored in forthcoming discussions. 

The spin-channel decomposition is presented in Fig. [3] for 
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FIG. 2. (a) Emca as function of number of Fe layers with (blue line) 
and without MgO (black line), (b) Emca as a function of Fe-0 dis¬ 
tance for 6-Fe/MgO. Atomic resolution of Emca versus number of 
Fe layer for the (c) clean Fe and (d) Fe/MgO 
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FIG. 3. (color online) Spin-channel decomposition of Emca for (a) 
clean-Fe and (b) Fe/MgO, where red (blue) symbols denote || (II) 
channel, and black ones for || channel. 


(a) the clean-Fe and (b) the Fe/MgO. Again, it is plotted as 
function of number of Fe layers, where the saturation behav¬ 
ior is evident. Remarkably, in both cases— with and without 
MgO, the || channel dominates over the other channels since 
the majority spin states are almost fully occupied as in the case 
of Fe multilayers 119, . 

To elucidate the physics origin of PMCA, band structures 
of Fe ML without and with MgO are shown in Fig. [4l We first 
analyze the 1-ML case for simplicity, to clarify the physics 
without loss of generality. The overall bands look so similar in 
two cases. The presence of MgO affects d z 2 and d XZJZ orbitals. 
Small shifts of d XZJZ are apparent due to MgO. At X , d xz are 
split into two states, one 0.6 eV below Ep and the other near 
Ef . On the contrary, in-plane orbitals are less affected by the 
presence of MgO. 

For the MCA analysis, we follow the recipe by Wang 
et al. Ill 911 . which has proved to be successful in various 
materials lf8L Fl9l425 1. The increase of MCA by the presence of 
MgO can be viewed in two ways. First, positive contribution 
by (xz\iz\yz) remains the same in spite of litte reduced contri¬ 
bution by (.x 2 —y 2 \£z\xy) around T. Second, negative contri¬ 
butions by (z 2 \£ x \yz) are reduced as the occupied d yz state be¬ 
comes unoccupied about \M-T. Third, negative contributions 
such as (xy\lx\xz,yz) in X-M are reduced owing to enlarged 
energy differences. As such, the presence of MgO slightly 
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FIG. 4. (color online) Band structure of the free-standing Fe (1-ML 
Fe) and Fe/MgO. Only the minority spin bands are shown, (a) 1-ML 
Fe,and (b) 1-ML Fe/MgO. The orbital contribution of Fe d state is 
emphasized in colors: red (d z 2), black {d x 2 _ y i), blue (d xy ), orange 
( d xz ), and green (d yz ), respectively. 


affects some bands due to weak hybridization. Nonetheless, 
the overall bands do not change much with a little increase in 
MCA. 

For completeness, the case of 6-ML Fe/MgO are analyzed. 
Contributions from Fe(I) and Fe(S) to bands are shown in 
Fig. 0a) and (b), respectively. Despite the complexity of 
bands, it is noticeable that the energy levels of those orbitals 
involving matrices (x 2 — y 2 \^z\xy) and (xy\£ x \xz,yz) do not 
differ very much in both Fe(I) and Fe(S). This implies that 
the hybridization is so weak that MCA from Fe(I) and Fe(S) 
are comparable. The large PMCA in Fe/MgO is a result of 
interplay of iz and t x matrices. 

To get more insights on the role of Fe(S) and Fe(I), DOS 
of the free-standing 6-ML Fe and the 6-ML Fe/MgO are pre¬ 
sented in Fig. 0c-f), where the left panel is for Fe(S) of the 
clean Fe and on MgO, and the right panel for Fe(I) and 0(1) in 
Fe/MgO. DOS of Fe(S) is little affected by MgO. The major¬ 
ity spin states of Fe(I) and Fe(S) are almost filled with peaks 
from d x 2_ y 2 and d z 2 1-2 eV above Ef. In particular, Fe(I) 
peaks from d z 2 in the majority spin state are closer to Ep. On 
the other hand, minority spin d^ and d XZJZ states close to Ep 
form peaks. The minority spin states of Fe(I) retain almost the 
same feature of the free-standing 6-ML Fe. DOS of 0(1) is 
also shown in Fig. 0f), where p states are prominent around 
-6 ^ -4 eV, which is rather far away from Fe d xzyz and d z 2 
states. As clearly seen from the DOS plots, the hybridization 
is weak in Fe/MgO. 

While MCA shows convergent behavior with the increase 
of the number of Fe layers, the singular behavior of 2-ML is 
puzzling. To tackle this issue, the unique structural feature of 
2-ML should be emphasized. Without MgO, both Fe layers 
are symmetrically equivalent. On the other hand, with MgO, 
while one layer is the surface layer the other is the interface 
layer. Moreover, the 2-ML exhibits rather unique electronic 
structure [See Supplementary Information]: d x 2_ y 2 state are 
prominent in DOS for both spin states just below Ep. These 
states couple with the minority-spin d XZJZ states, hence pos¬ 
itive contribution is compensated by negative one. Also, Fe 
2pi/ 2 core-levels are analyzed. The core-level shifts reflect 
either charge transfer or change of internal electric field i27ll. 
Regardless of the presence of MgO, 2 p l / 2 energy of the center 
Fe layer is lowest, while those of the interface and the surface 
are higher. [See Supplementary Information for core-levels.] 
In other cases, Fe layers are thick enough to screen out in¬ 
ternal field by MgO, which is not the case for 2-ML Fe, as 
manifested in Fe 2 pi/ 2 core levels. 

So far in our calculations, the stoichiometry is assumed to 
be perfect. On the other hand, in over- and under-oxidized 
Fe/MgO[9|], the Fe layer in contact with the over- or under¬ 
oxidized MgO is expected to show ripples Here, we argue 
that perfect epitaxy at the interface is the key to achieve large 
PMCA in Fe/MgO. Results clearly show that MCA contribu¬ 
tions from Fe(S) and Fe(I) are similar in magnitudes regard¬ 
less of the hybridization. Hence, the hybridization is weak as 
evidenced in DOS, which cannot be the main driving force for 
large PMCA. although its contribution to PMCA cannot com- 
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(a) Fe (I) in Fe/MgO (b) Fe (S) in Fe/MgO 
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FIG. 5. (color online) Band structure of (a) Fe (I) and (b) Fe (S) in 
6-ML Fe/MgO. Color representing d orbitals are the same as Fig.|4] 
Density of states (PDOS) of Fe and Fe/MgO with 6-ML Fe. The 
surface Fe (c) without MgO, (d) Fe(I), (e) Fe(S)with MgO, and (f) 
0(1) 2 p orbitals. The d orbital states are shown in different colors: 
red (d z 2), black (d x 2_ y i), blue ( d xy ), and green (d XZJZ ), respectively. 
For the 2 p orbitals: p x , Py, and p z , respectively. 


pletely neglected. Before we conclude, we emphasize here 
that the physics we have demonstrated occurs in other lattice 
constants as well. Even if we cannot exclude the strain effect 
in Fe/MgO, we argue that the main driving force of PMC A is 
the perfect epitaxy which is intrinsically achieved in ab initio 
calculations. 

In summary, we have shown that the p-d hybridization, 
though not completely ignorable, cannot be the main driv¬ 
ing force for large PMCA, which is different from commonly 
accepted. The hybridization is overall weak, and indeed in¬ 
creases systematically PMCA with MgO. However, MCA 
with and without hybridization are almost equal as well man¬ 
ifested in Fe(S) and Fe(I). From the detailed analysis, contri¬ 
butions from the surface and the interface Fe are dominant. 
As in the Fe multilayers, since d states of the majority spin 
bands are completely filled, it is the | j. channel that plays the 
most dominant role in PMCA. Furthermore, while the pres¬ 
ence of MgO retains positive contribution by (xz\£z\yz), neg¬ 
ative contribution (z 2 \£x\yz) and (xy\£x\xz,yz) are slightly re¬ 
duced. More importantly, perfect epitaxy, which is an indica¬ 
tor of high level of crystal growth, is the key factor to realize 
large PMCA in Fe/MgO from the fact that Fe layer with and 


without hybridization contribute almost equally to PMCA. 
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